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PREFACE 


The  present  Research  Memorandum  presents  a  general  computer  program 
developed  In  the  Statistics  Research  and  Analysis  Laboratory  of  the  U.  S. 
Army  Personnel  Research  Office  for  use  in  the  analysis  of  simulated  man¬ 
power  information  systems  problems  having  common  aspects. 

With  feasible  modifications ,  the  program  is  applicable  to  a  number 
of  major  studies  of  manpower  allocation  policy  and  procedures. 
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A  GENERAL  COMPUTER  SIMULATION  FCR  CONDUCTING  ALLOCATION  EXPERIMENTS 


L..  studies  of  tie  Arab's  pereacnel  systems,  t'.af  effects  ;  -posed 
alternative  policies  and.  procedures  must  be  evslvT-d.  m  lie  a  or  ex¬ 
perimentation  vith  the  real  system- -trial  isspleac.r&ation  of  proposed 
changes  ,  for  example — the  U.  S.  Army  Personnel  Research  Office  haa  de¬ 
veloped  a  series  of  models  representing  the  system.  Some  of  these  models 
have  been  Computerized  and  made  dynamic — computer  simulations,  in  other 
words . 

A  number  of  experiments  have  been  performed  in  U.  S.  AHlO’s  Sta¬ 
tistical  Research  and  Analysis  Laboratory  in  which  simulated  samples  of 
enlisted  men  were  generated  and  treated  experimentally.  Optimal  alloca¬ 
tion  to  job  clusters  associated  vith  the  eight  Army  aptitude  areas  on  the 
basis  of  performance  estimates  linearly  transformed  from  scores  on  the 
Army  Classification  Battery^*' was  an  important  feature  of  the  experi¬ 
ments.  In  one  such  study,  the  transformation  consisted  of  least  squares 
regression  coefficients  based  on  all  eleven  tests  of  the  ACB  rather  than 
on  computationally  simplified  two-test  composites — the  aptitude  area 
scores.  It  was  found  that  the  sum  of  the  performance  estimates  for  the 
jobs  to  which  the  men  were  assigned  was  substantially  increased  by  the 
experimental  procedure.  This  "allocation  sum"  has  also  been  studied  as 
a  function  of  the  metric  in  which  the  performance  estimates  were 
expressed^.  For  example,  performance  estimates  based  on  100-unit  equal 
interval  scales  were  found  to  yield  higher  overall  estimated  performance 
than  estimates  based  on  ordinal  scales  constructed  by  ranking  performance 
on  each  job  over  individuals  or  by  ranking  each  individual  over  Jobs. 

Each  of  these  scales  was  found  to  be  superior,  however,  to  equal  interval 
scales  vith  only  10  units.  In  studies  with  a  still  different  orientation, 
predicted  changes  in  overall,  performance  were  examined  when  different 
screening  or  allocation,  policies  for  enlisted  men  were  simulated.  Alter¬ 
native  policies  consisted  of  variations  in  requisites  for  assignment  to 
occupational  areas  in  terms  of  estimated  performance  levels.  Estimates 
were  based  on  scores  on  the  Armed  Forces  Qualification  Test  (AFQT),  the 
Army  Qualification  Battery  (AQB),  the  Army  Classification  Battery  (ACB), 
or  sane  specified  combination  of  scores. 


^Sorenson,  Richard  C.  Optimal  allocation  of  enlisted  men — Full  regres¬ 
sion  equations  vs  aptitude  area  scores.  USAPRO  Technical  Research 
Note  165,  November  1965. 

^Sorenson,  Richard  C.  Effect  of  modification  in  operational  procedures 
on  an  optimal  personnel  allocation  system.  Paper  presented  at  the 
Convention  of  the  American  Psychological  Association,  September  1965. 
Abstract  appeared  in  American  Psychologist  20;7,  1965. 

^Sorenson,  Richard  C.  This  effect  of  metric,  changes  on  resource  alloca¬ 
tion  decisions.  Proceedings  of  the  U.  S.  Army  Operations  Research 
Symposium,  Fort  Monmouth,  New  Jersey,  March  1966. 


Each,  of  these  ex, percents  has  required,  the  preparation  of  a  separate 
and  lengmiy  computer  ^  vrant„  Yet,  there  was  cor-siderable  similarity  in 
the  programs .  To  facilitate  experiments  involving  the  generation  of 
simulated  samples,  U.  S.  1PP0  research  scientists  have  developed  a  gen¬ 
eral  simulation  program  which  may  be  used  to  evaluate  a  varieV  of  alter¬ 
native  policies.  The  program,  designed  to  perform  the  many  operations 
common  to  experiments  of  this  kind,  can  easily  be  modified  to  handle 
computations  specific  to  different  experiments. 


general  operations 

The  general  part  of  the  program  begins  with  the  automatic  generation 
of  a  vector  of  random  normal  deviates  for  every  individual]  on  this  vector 
are  performed  a  series  of  linear  transformations  of  the  form  V  =  JJ  T  +  M. 

T  may  be  a  matrix  of  least  squares  regression  coefficients  for  obtaining 
performance  estimates  V  from  a  set.  of  predictor  scores  U.  M  are  the 
additive  constants  to  yield  estimates  with  specified  means.  A  special 
purpose  for  which  T  is  used  it  simulation  studies  is  to  transform  randcm 
normal,  deviates  which  have  an  expected  covariance  matrix  of  identity 
into  variates  with  an  expected  covariance  matrix  characteristic  of  the 
population  under  investigation.  The  user  specifies  the  series  of  linear 
transformations  required  to  generate  a  particular  sample  by  first 
inputting  the  matrices  (and  associated  mean  vectors)  and  then  referencing 
these  matrices  on  special  transformation  cards,  which  are  input  in  the 
order  of  the  transf ormations  tc  be  perforated 


SPECIALIZED  SUBROUTINES 


To  perform  non-linear  transformations  on  the  score  vectors  generated 
for  each  individual,  specially  written  parameter  subroutines  may  be  easily 
incorporated  in  the  sequence  of  operations  under  control  of  the  main  pro¬ 
gram.  Modification  and  recompilation  is  required  only  for  a  short  execu¬ 
tive  subprogram,  net  for  the  main  program.  The.  parameter  subroutines  are 
assigned  integer  names  and  are  called  by  listing  these  integers,  in  the 
order  in  which  the  subroutines  are  to  be  performed,  on  the  cards  which 
define  the  sequence  of  linear  transformations.  These  subroutines  are  also 
used  to  determine  whether  the  scores  being  generated  for  a  given  individual 
axe  consistent  with  sample  characteristics  defined  for  a  particular  in¬ 
vestigation.  As  the  operations  specified  on  each  transformation  card  are 
completed,  a  special  parameter,  which  may  have  been  set  to  Indicate  re¬ 
jection  of  the  given  Individual  by  any  of  the  parameter  subroutines,  is 
automatically  sensed  by  the  main  program.  Thus,  tests  for  individual 
acceptance  or  rejection  may  be  performed  repeatedly  and  at  any  stage  in 
the  computations. 

The  result  of  the  sequence  of  operations  specified  on  the  transforma¬ 
tion  cards  is  the  construction  of  a  v  x  p  matrix  of  performance  measures) 
the  P  columns  correspond  to  Job  categories  under  investigation,  and  the  v 
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rows  represent  the  individuals  in  the  sample.  Based  on  this  matrix, 
each  of  the  v  individuals  is  assigned  to  one  of  the  P  jobs  in  such  a  way 
that  required  quotas  far  the  different  jobs  are  filled,  and  the  sum  of  the 
performance  estimates  for  the  jobs  assigned  to  the  v  individuals  is 
maximal.  The  solution  for  this  optimal  allocation  may  be  consxoered  a 
special  case  of  Kuhn’s4'  solution  to  tne  transportation  problem. 

After  optimal  allocation,  the  response  vectors  for  the  v  individuals 
are  regenerated  to  compute  statistics  that  are  functions  of  the  job  to 
which  each  individual  is  assigned.  The  programmer  specifies  the  new  com¬ 
putations  by  inputting  a  second  set  of  transformation  cards  of  the  type 
input  before  allocation.  (To  eliminate  unnecessary  computations,  special 
codes  on  the  transformation  cards  can  be  used  to  restrict  arithmetic 
operations  to  those  involving  the  assigned  jobs.)  The  same  sequence  of 
random  numbers  is  automatically  generated  both  before  and  after  allocation 
so  that  the  simulated  performance  measures  will  represent  the  same 
individuals . 

The  feature  of  specifying  alternative  sequences  of  linear  transforma¬ 
tions  and  subroutines  as  parameters  may  be  considered  a  compiler  for  in¬ 
dividual  response  simulation.  Unfortunately,  incorrect  use  of  this  part 
of  the  program  will  not  produce  simple  error  diagnostics,  but  may  generate 
meaningless  results.  To  insure  that  the  same  individuals  are  accepted 
into  the  sample  before  and  after  allocation,  the  user  must  specify  the 
correct  sequence  of  transformations  and  subroutine  calls.  In  addition, 
care  must  be  taken  to  maintain  the  correct  order  of  simulated  scores  which 
represent  tests  of  interest.  For  example,  the  transformation  V  =  U  T, 
which  can  be  specified  on  a  transformation  card  by  equating  an  integer 
to  the  order  in  which  the  appropriate  matrix  was  input,  will  operate  only 
on  the  first  a  elements  in  the  first  working  storage  area  U  and  wilj^  re¬ 
place  only  the  first  0  elements  of  the  second  working  storage  area  V; 
a  and  0  are  the  number  of  rows  and  columns,  respectively,  of  the  trans¬ 
formation  matrix  T.  Scores,  then,  which  are  to  be  modified  by  reference 
to  parameter  subroutines,  but  not  by  post -multiplication  by  one  of  the 
transformation  matrices,  must  be  represented  by  elements  subsequent  to 
those  altered  by  linear  transformations.  Accuracy  of  results  can  be 
checked  by  setting  a  special  debug  parameter,  which  will  produce  detailed 
print-out  of  computations  intermediate  to  the  final  results. 


THE  COMPUTATION  PROCEDURE 

The  program  operates  in  four  parts :  I,  input  and  initialization  re¬ 
quired  before  sample  generation  II,  simulation  of  individual  response 
vectors,  optimal  allocation,  and  regeneration  of  response  vectors;  III, 
computations  and  output  based  on  the  v  individuals  in  a  sample,  plus 


^  Kuhn,  H.  W.  The  Hungarian  method  for  the  assignment  problem.  Naval 
Research  Logistic  Quarterly,  2:1,  1955« 
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initialization  required  for  generation  of  a  new  sample;  and  IV,  computa¬ 
tions  and  output  based  on  all  samples  generated-  Reference  to  the  flow 
diagram,  Figure  1,  may  clarify  this  sequence  of  operations. 


To  perform  operations  consistent  with  the  four  parts  of  the  main 
program,  four  entry  points  must  be  written  into  every  parameter  subroutine. 
To  transfer  control  to  parts  I,  III,  and  IV,  the  main  program  automatically 
calls  all  subroutines  listed  on  the  transformation  cards  in  order  of  in¬ 
creasing  size  of  integer  names  assigned  to  the  subroutines,  (During  execu¬ 
tion  of  part  II,  the  subroutines  are  called  in  the  order  listed  on  the 
transformation  cards.)  Dummy  instructions  ("returns”)  are  written  into 
entries  for  any  of  the  four  parts  net  used  by  a  subroutine. 


INPUT  OPERATIONS 

Operations  which  constitute  Part  I  consist  of  the  input  of  trans¬ 
formation  matrices  and  associated  mean  vectors  (box  l),  input  of  the 
vector  of  starting  numbers  to  initialize  random  number  generation  (box  2), 
input  of  the  number  of  individuals  to  be  simulated  and  other  parameters 
defining  characteristics  of  each  sample  (besx  3 ) >  input  of  ,ne  transforma¬ 
tion  cards  (box  6)  and  a,  the  number  of  samples  to  be  simulated  (box  7). 

To  make  it  convenient  to  generate  6  different  sets  rf  samples  during  one 
computer  run,  each  set  being  defined  by  different  lists  of  transformation 
cards,  the  parameter  5  is  input  (box  4).  All  transformation  matrices  re¬ 
quired  for  the  6  sets  of  samples  are  input  simultaneously  under  control 
of  the  first  part  of  the  program  (box  l).  Any  of  these  matrices  may  he 
‘used  by  any  of  the  6  different  problems  by  entering  an  integer  which 
corresponds  to  the  order  in  which  the  given  matrix  was  input  into  the 
appropriate  column  of  a  transformation  card. 

The  input  of  parameters  of  data  required  for  execution  of  any  of  the 
parameter  subroutines  is  represented  by  box  8.  This  input  is  called  in 
order  of  increasing  size  of  the  integer -named  subroutines.  Any  additional 
initialization  required  for  subsequent  operation  of  the  subroutines  is 
also  performed  within  this  first  entry. 


SIMUIATION  OPERATIONS 

Part  II  of  the  program  is  illustrated  in  boxes  9-27  of  the  flow 
diagram.  In  boot  the  starting  vector  Afln#  which  controls  the  generation 

of  random  numbers,  is  set  equal  to  the  vector  input  to  represent  the  first 
individual  of  sample  1  (A^.  in  box  2).  This  vector  automatically  changes 

as  random  numbers  are  generated  for  each  individual  (boxes  21  and  27). 

Bax  16  refers  to  the  case  where  the  responses  being  generated  for  indi¬ 
vidual  i  have  been  found  unacceptable  by  one  of  the  parameter  sub¬ 
routines;  consequently,  another  vector,  A  ,  is  generated  from' which  a 

“'Dll 
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Figure  1.  Flow  diagram  of  General  Program  for  Conducting  Optimal  Allocation  Experiments 

PART  I  -  Input  Operations 


Figure  1.  (continued)  PART  II  -  Simulation  Operations 
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Figure  1  (continued). 

PART  III  -  Computations  and  Output  (  v  individuals  )  and 
PART  IV  -  Computations  and  Output  (  a  samples) 
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new  set  of  scores  representing  the  n  individual  can  be  constructed. 

To  insure  that  the  same  set  of  individual  responses  are  generated  before 
and  after  allocation,  t.r.e  starting  vector  for  individual  1  of  sample  s, 


A  ,, 
-si* 


is  saved  for  the  re -initialization  shown  in  box  23  •  The  starting 


vector  which  defines  the  beginning  of  each  new  sample  is  autcnsar.i- ally- 
printed  out  by  the  main  program  (box  10)  so  that  simulation  can.  be 
continued  from  any  point  of  unanticipated  termination  of  production. 


Linear  transformations  and  parameter  subroutines  specified  by  the 
first  t  transformation  cards  are  performed  on  each  individual  response 

vector  of  random  normal  deviates  (boxes  13  and  14)  until  the  matrix  cor 
sisting  of  c  performance  measures  for  v  individuals  has  been  constructed. 
Execution  of  the  transformation  cards  is  under  control  of  the 

index  t,  as  shown  in  boxes  17-19,  21,  23-25,  27.  Box  15  illustrates  the 
feature  of  examining  the  qualifications  of  each  individual  every  time  the 
computations  specified  on  a  transf ormation  card  are  completed. 


COMPUTATIONS  AND  OUTPUT 

Parts  III  (box  28)  and  IV  (box  31)  represent  the  automatic  calls  to 
the  parameter  subroutines  to  perform  operations  over  the  v  individuals  in 
a  sample  and  then  over  the  a  samples  respectively.  The  final  loop  & 
determines  whether  a  new  set  of  transformation  cards  will  be  input  (at 
box  6)  to  define  additional  samples.  The  new  samples  will  be  based  on 
the  same  set  of  parameters  (v,  tt,  p,  boxes  1-4)  which  defined  the  previous 
samples,  with  differences  occurring  only  as  specified  on  the  new  trans¬ 
formation  cards.  Note  that  the  starting  vector  which  begins  the  genera¬ 
tion  of  random  numbers  is  reset  to  A^  (box  9),  so  that  equivalent  sets 

of  individual  responses  will  be  generated  for  each  of  the  6  sets  of 
samples . 


MANUAL  FOR  USE  OF  GENERAL  COMPUTER  PROGRAM 
FOR  SIMULATION  EXPERIMENTS 


INPUT  FORMAT 

iter  1  i  J  (Format*  1XA6,I4, 12) 

i  is  an  integer  which  occupies  columns  8-11  and  specifies  t.be 
number  of  complete  sets  of  data  and  parameters  to  be  input. 

To  debug  new  subroutines,  detailed  print -out  of  individual 
score  vectors  may  be  obtained  by  setting  i  =  1,  columns  12-13; 
otherwise,  j  =  0. 
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TRANSFORMATION  MATRICES 


lj,  2  j , . . . ,  i  j , . . ,  (For?..*-,  10(11,11, Al)) 

Beginning  in  column  1.  a  two-digit  integer  ij,  which  nay  be 
followed  by  a  conroa,  is  listed  for  every  matrix  input,  i  =  1 
for  the  first  matrix  input  and  is  increased  by  one  for  each 
successive  matrix;  the  maximum  i  =  nmat  ("nmat"  is  a  para¬ 
meter  set  within  the  main  program  to  determine  the  total 
number  of  transformation  matrices  to  be  input).  If  a  vector 
of  means  associated  with  a  given  transformation  is  to  be 
input,  j  =  1;  otherwise,  i  =  0. 

or  3  (Format  212) 

MATRIX  X  (Format  IX,  9A8) 

Transformation  Matrix  (Fermat  9^8 . 1 ) 

For  every  ij  listed,  one  card  is  input  defining  the  order  of 
the  matrix,  followed  by  a  second  card  containing  any  alpha¬ 
numeric  name  for  the  matrix,  followed  by  the  matrix  itself. 
a  (columns  1-2)  is  the  number  of  rows,  and  3  (columns  3“4) 
is  the  number  of  columns.  Elements  of  the  matrix  are  input 
row  by  row  with  the  format  9F8.4. 

* 

MATRIX  XM  (Format  IX,  9A8) 

Means  Following  the  transformation  matrix  associated  with  a  given 

ij,  if  j  =  1,  an  additional  alphanumeric  card  identifies  the 
vector  of  means.  The  means  begin  in  the  next  card  with 
format  9F8.4. 

Starting  Vector 

The  y  =  24  octal  digits  which  begin  the  generation  of  y 
independent  sequences  of  random  numbers  are  input  on  3  cards 
with  format  908  for  the  C.D.C.  computer  (and  on  4  cards  with 
format  6012  for  the  IBM  7094s-'  at  NBS). 


Canmercial  designations  are  given  in  the  interest  of  precision  in 
specifying  the  equipment  used,  and  do  not  constitute  official 
indorsement  by  the  Army.  * 


v  it  p  (Format  14,212) 

v  (columns  .1-!* )  ts  the  number  of  individuals  in  the  sample . 

tt  (columns  5-6) is  tbs  number  of  independent  randcec.  ccrmal 
deviates  generated  for  each  iadividual  to  initiate  toe 
simulation  of  p  performance  measures. 

p  (columns  7-8)  is  the  number  of  job  categories  on  which 
optimal,  allocation  is  to  be  based. 


J*1  (yowaat  llfI5) 


Required  quotas  for  each  of  P  jobs:  £  «  v. 

Jal  J 


iter  2  -  6  (Format  1XA6,I4) 


ft  specifies  the  number  of  complete  sets  of  parameter  cards 
and  data  which  are  to  follow,  all  of  which  will  be  controlled 
t»y  the  set  of  cards  previously  input. 


TRAflSFCRNATION  CARDS 
T  T2  (Format  212) 

cards  will  follow  specifying  transformations  and  sub¬ 
routines  to  be  performed  on  the  vector  of  n  normally  dis¬ 
tributed  standard  scores  generated  for  individual  n  at 
working  storage  U,  preceding  optimal  allocation.  The  next 
t2  cards  specify  the  transformations  and  subroutines  to  be 

performed  on  the  same  individual  response  vector  regenerated 
at  U  following  optimal  allocation. 

The  post -multiplication  of  a  response  vector  located  at 
working  storage  U  (or  V)  by  any  transformation  matrix  Input 
(Tj),  the  addition  of  the  vector  of  means  associated  with  the 

given  transformation  (Mj)>  the  storage  of  the  result  in 

working  storage  V  (or  U)  can  be  specified  on  any  of  the 
+  t2  transformation  cards  by  the  following  codes: 


12,1  =$► 

V 

=  U 

^1 

+  M  j 

lx?! 

lXOj 

«IX 

h 

3xP. 

21,1  => 

U 

=  V 

+  M  j 

lxgj 

aIx 

lap. 
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Different  Tj  and  are  referenced  by  setting  I  to  any  integer 

from  1  to  nmat  corresponding  to  the  order  in  which  the  trans¬ 
formation  matrices  are  input.  If  no  vector  of  means  has  been 
imput  for  1T,  ML  =  0.  Multiplications  are  performed  on  only 

the  first  or.  elements  of  U  (or  V),  and  the  result  changes  only 

the  first  pj  elements  of  V  (or  U);  ctj  and  are  the  number  of 

rows  and  columns  of  T^. 

Y  scores  can  be  relocated  from  one  working  storage  area 
to  the  other  by  setting  1=0.  That  is. 


\  ew 

12,0  =$>  V  =  U 

lxY  IxY 


21,0  =#►  U  =  V 
IxY  IxY 


Following  allocation  over  jobs  j  =  1,  2,  .  .  .,  p,  a 
statistic,  which  is  a  function  of  the  Job  J  to  which  egch 
individual  is  assigned,  may  be  computed  and  stored  in  x 
(W(k9)).  That  is. 


10,1 


x 

lxl 


U 

lxa. 


T 

-Ij 

«jX  1 


+  m. 


IJ 

lxl 


20,1 


x 

lxl 


V 
lxQf 


T 

I  OfjXl 


m. 


Ij 

lxl 


where  is  the  j-th  column  of  and  m^j  is  the  j-th 
element  of  M_.  If  no  transformation  matrix  is  to  be  used, 

rs*  *** 

the  j-th  element  of  U  (or  V)  can  be  relocated  in  x  by 
setting  1=0: 
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PARAMETER  SUBROUTINES 

To  perform  additional  kinds  of  computations  on  scores  in 

/■w  /V 

U,  V,  or  x,  a  list  of  subroutines  may  follow  the  transforma¬ 
tion  codes  on  any  of  the  t^  +  t^  transformation  cards .  Each 

subroutine  is  equated  to  some  'integer  between  1*  and  99  and  is 
listed  in  the  ^rder  of. required  execution.  ’Beginning  in 
column  5  with  format  5(A1,I2),  each,  subrouting  is  introduced 
by  a  comma,  with  the  following  integer  n£me  occupying  two 
columns.  As  many  as  5  subroutines  may  be  referenced  on  a  • 
given  card.  To  illustrate,  suppose  the  post-multiplication 

.f of  scores  in  U  by  the  second  matrix  inpvrc  is  required,  to  be 
■followed  by  the  execution  of  subroutines*  3/ 'S&>  and  1,  which 
have  been  especially  written* to  Operate  on  the  result  in 

ew  I 

V.  The  transformation  card  would  appear' as: 

12,2,03,98,01 

During  simulation  of  scores  for  a  giVen  individual,  a 
reject  parameter  is  evaluated  each  time  the  processing  de¬ 
fined  by  a  transformation  card  is  completed.  Depending  on 
the  value  of  this  parameter,  .scores  for  the  individual  are 
either  retained  for  further  processing .or  are  treated  as  in¬ 
appropriate  for  a  given  sample  and  eliminated.  If  it  is 
desired  to  subdivide  a  list  of  subroutines  so  that  the  re-  ' 
ject  parameter  will  be  evaluated  after  the  execution  of  only 
one  or  two,  say,  of  the  subroutines,  additional  transforma¬ 
tion  cards  defining  the  subroutine  subsets  are  inserted,  but 
with  columns  1,  2,  and  4  set  to  zero.  To  illustrate,  suppose 
the  y  scores  in  V  are  to  be  transferred  to  U,  to  be  followed 
by  execution  of  subroutines  3,  98,  and  1,  written  to  operate 
on  scores  in  U.  The  reject  parameter,  modified  by  subroutines 
3  and  1,  say,  is  to  be  evaluated  after  the  execution  of  3  and 
1.  The  appropriate  two  transformation  cards  would  appear  as: 

21,0,03 

00,0,98,1 


NUMBER  CF  SAMPIES 

iter  1  =  i  (Format  1XA6,I4) 

i  is  an  integer  which  occupies  columns  8-11  and  specifies  the 
number  of  independent  samples  to  bs  simulated  as  defined  by 
all  previous  parameter  cards  input. 
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INPUT  CALLED  BY  SURBOUTINES 


Tbs  input  ex'  parameters  and  data  may  be  required  for  the 
execution  cf  subroutines  listed  on  the  transformation 
cards.  Input  called  by  different  subroutines  is  arranged 
by  increasing  order  of  the  integers  which  name  the  sub¬ 
routines  . 


DEFINITION  CF  DIMENSION  STATEMENTS  AS  LISTED 
IN  THE  COMPUTER  PROGRAM 


KO(N,LC)s  Mat  rix  of  fLced-pcmt  performance.  estimates  for  N  individuals 

on  each  of  1C  jets.  Currently  (February  1967 ),N  ~  1000  and  LC  =  8 

IROW(N):  Vector  containing  the  job  j ,,  j  -  1,  2,  .  .  .  ,  LC,  to  which 

each  of  N  individuals  is  assigned. 

Kftl(LC),  KQ2(LC): 

Storage  areas  for  the  required  quotas  for  each  of  LC  Jobs. 

Z(AREA):  Storage  area  fer  all  transformation  matrices,  thus  requiring 

that 

NMAT 

AREA  =  L  (NA  x  MA  ), 
m=l  31 

where  the  order  of  each  matrix,  input  is  NA  (rows)  x  MA 
(columns),  and  the  total  number  of  matrices  input  is  NMAT. 

Currently,  AREA  =  (23  x  23)  +  (11  x  8)  +  (11  x  8)  -  703. 

ZM(TMA) :  Storage  area  for  means  associated  with  each 'transformation 

>  NMAT 

matrix,  with  TMA  =  L  MA 

m=l  m 

Currently,  TMA  =  23  +  8  +  8  =  39. 

U(S)^,  V(S),W(2S  +  1),  with  U(l)  =  W(l)  and  V(l)  -  W(S  +  l): 


The  storage  area  represented  bare  by  S  is  represented  in  the  descrip* 
tion  of  the  program  by  Y« 

-  13  - 


1 


f  Working  storage  areas  a? ex  for  siEuletirv  a  maximum  of  S 

i  *'.»  ^  i*J\  l  1  :•  ;:  : 

[  *•  -•>«.«•  1  uuvk  -«..A  &  —  ..  o 

t  S' 

|  Ourrentu.y  .February  lye  ~ ,  S  •=>  >k. 

IISV(S,\  ISV(S),  JSV(S)j 

Areas  use.*,  for  storf  .*  the  starting  vr  cor  of  S  octal  digits 
I  called  by  -.ufcrouf  ine  RAKD. 

BBB(S)  (required  for  the  ISM  7094': 

Working  storage  area  used  by  sub roue in-  RAND. 

VF0R(9)  (for  CDS)  or  WCR{l£)  (for  IBM. 

Area  used  for  input  of  variable  format. 

IS(NTRAN),  JS(HSRA»,5).»  •IN(NTPAM)J,  K2f(JERAN,6).  L35l(JSS»JB)i 

Storage  areas  used  for  setting  up  sequence  of  transform  - 
;  tions  and  subroutines  used  in  simulating  responses  for  each 

?  individual.  RERAN  1?  the  maximum  number  of  transformation 

;  •  cards  which  can  be  input  for  a  given,  simulation  experiment. 

NSUB  is  tbs  maximum  number  of  different  parameter  subroutines 
which  can  be  listed  on  trass format ion  cards  for  tie  experiment. 
;  Currently,  NIRAIJ  =  8  and  NSUB  =  10. 

;  na(nmat),  ma(nmat),  rb(nmat } : 

j 

•  Areas  used  for  addressing  brans  formation  matrices  and 

associated  vectors  of  means.  Currently,  NMAT  *  J. 


For  tbs  most  efficient  utilization  cf  ecaaputer  space,  tie 
main,  program  can  be  recompiled  so  that  N,  LCt,  AREA,  TMA,  and 
NMAT  are  consistent  with  requirements  for  specific  projects. 

The  statement  "NMA1  =  3,”  which,  is  within  the  first  part  of 
tbe  main  program,  will  also  have  to  be  changed.  Changes  in  S, 
RERAN,  and  NS’JB  will  not  result  ir.  appreciable  saving  cf  space. 
If  change?  are  required  for  S,  NTRAI7,  and  NSUB,  additional 
changes  must  be  made  in  statements  of  the  main  program. 


PROGRAMMIN';  RECUIREMBWTS  FOR  PARAMETER  fFlRFOUTTNES 


Any  subroutine  wr :  ~>:..er.  to  ha.udlr  computations  specific  to  different 
projects  must  be  assoc  rated  w.tn  an  integer  .namr,  and  the  first  parameter 
of  every  subroutine  must  be  "JSUB".  Thus,.  "SUBROUTINE  SUB  J  (JSJB,  ...)", 
where  J  is  any  integer,  would  appropriately  been-,  a  subroutine.  In  addi¬ 
tion,  the  first  statement  of  any  subroutine  must  be  "GO  TO  (a,  b,  c_,  d), 
JSUB",  a,  b,  c,  and  d  not  all  necessarily  different. 

All  subroutines  listed  on  the  transformation  cards  are  automatically 
called  in  four  different  parts  cf  the  main  program.  Corresponding  to 
each  part,  .JSUB  is  set.  to  1.  2,  J.f  or  4  to  permit  execution  of  different 
types  of  subroutine  operations. 

Before  beginning  response  Simula  lor  ary  sample  of  w  individuals, 
the  main  program  sets  JSUB  *  I  and  transfers  : octroi  to  statement  a  of 
each  subroutine  in  order  of  ::n.:r.  easing  size  of  the  integers  naming  the 
subroutines.  Ary  initialization  or  calls  for  input  of  data  or  parameters 
required  for  subsequent  operations  must  be  written  into  this  "part,  a"  of 
each  subroutine.  If  no  such  initialization  is  required,  statement  a  will 
be  a  "RETURN". 

During  the  simulation  cf  individual  score  vectors  in  the  second  part 
of  the  main  program,  JSUB  =  2;  for  earn  individual  in  every  sample,  con¬ 
trol  passes  to  statement  t  of  each  subroutine  according  tc  the  order  trans¬ 
formations  and  subroutines  are.  listed  on  the  transformation  cards.  The 
main  function  of  subroutine  sections  beginning  at  statement  b,  then,  will 
he  to  operate  on  scores  in  wording  storage  areas  tf,  V,  or  x  or  to  construct 
KO,  the  v  x  p  matrix  of  scores  used  by  subroutine  OTT  for  optimal  alloca¬ 
tion. 

After  all  computing  to  be  performed  on  v  individuals  has  been  com¬ 
pleted,  the  main  program  sets  JSUB  =  j,  and  control  is  transferred  to 
statement  c_  of  each  subroutine,  again  in  order  of  Increasing  size  of  the 
integers  naming  the  subroutine.  "Part  c"  of  each  subroutine  must  be 
written  to  perform  any  summary  computations  based  on  all  v  individuals 
for  the  sample  just  completed,  and  to  perform,  any  initialization  required 
for  generating  the  next  sample.  If  no  such  computations  are  required, 
statement  c  must  be  a  "RETURN" . 

When  simulation  of  the  total  of  a  samples  has  been  completed,  the 
main  program  sets  JSUB  =  4,  and  control  is  transferred  to  statement  d  of 
each  subroutine,  again  according  to  the  order  of  the  subroutine  numbers. 

Any  summary  computations  or  output  based  on  the  completed  set  of  o  samples 
or  any  initialization  required  before  parameters  defining  a  new  set  of  a 
samples  dre  input  must  be  written  beginning  at  statement  d.  Otherwise, 
statement  d  must  be  a  "RETURN". 
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MODIFICATION  OF  SUBROUTINE  EXSUB 


The  call  statement  for  any  of  the  parameter  subroutines  listed  on 
the  transformation  cards  occurs  in  the  executive  subroutine.  Each  time, 
then,  that  a  new  subroutine  is  prepared  for  a  specific  project,  "Sub¬ 
routine  EXSUB"  must  be  recompiled. 

At  the  time  for  execution  of  any  of  the  subroutines  listed  on  the 
transformation  cards,  the  main  program  automatically  sets  a  parameter 
"KM"  within  "EXSUB"  to  the  integer  name  of  the  subroutine.  Any  test 
which  recognizes  that  KM  equals  this  integer  name,  J,  say,  and  is 
followed  by  the  execution  of  a  statement  calling  "SUBROUTINE  J",  may  be 
written  into  the  executive  routine.  One  convenient  method  would  be  to 
extend  the  length  of  the  computed  go-to  statement  (currently,  "GO  TO 
(1,2, ?, 4, 5), KM"),  making  the  J-th  parenthesized  integer  equal  to  J  such 
that,  when  KM  =  J,  control  will  be  transferred  to  statement  J,  and  sub¬ 
routine  J  will  be  executed.  Care  must  be  taken  to  insert  any  new  sub¬ 
routine  call  statements  to  operate  within  the  range  of  the  do-lopp  which 
terminates  at  statement  (currently  labeled)  25. 

All  regions  likely  to  be  referenced  by  specific  subroutines  are  in 
the  COMMON  area  for  both  the  main  program  and  EXSUB.  Those  regions  most 
likely  to  be  needed  as  parameters  in  any  of  the  subroutines  are  defined 
as  follows : 

KO  the  v  x  p  matrix  of  fixed-point  scores  used  by  subroutine 

OTT  for  optimal  allocation. 

U  and  V  working  storage  areas  beginning  at  W(l)  and  W(2?),  respectively, 
used  for  the  post-multiplication  of  an  individual  score  vector 
by  a  transformation  matrix  and  possible  subsequent  addition  of 
a  mean  vector. 

LC  (equivalent  to  p) 

a  single  score  selected  from  the  p  scores  simulated  for  a 
given  individual  to  correspond  to  j,  j  ■  1,2,..., p,  the  job 
to  which  the  individual  is  assigned. 

N  the  number  of  Individuals  in  a  given  sample. 

I  index  for  the  number  of  individuals ;  I  ■  1,2, ...,N 

•• 

ITER?  (equivalent  to  3) 

the  number  of  samples  of  "N"  individuals  to  be  simulated  for 
each  set  of  transformation  cards  input. 

IT?  index  for  the  number  of  samples:  IT?  »  1,2,..., ITER? 
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